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Abstract The fundamental part of a dowel-type connec¬ 
tion is the embedment of a steel dowel into the timber that 
surrounds it, and the stiffness of the timber in embedment 
is represented by the foundation modulus. A standard 
experimental method for identifying the foundation mod¬ 
ulus under static load is modified to assess the secant 
stiffness exhibited under one-sided cyclic load. It is shown 
that the steady-state secant stiffness is significantly higher 
than the static stiffness under initial loading, and that, if the 
amplitude of the cyclic component of the load is suffi¬ 
ciently small, a simple analytical elastic model predicts the 
foundation modulus well. The analytical model is based on 
a complex stress function for the timber in embedment and 
the frictional interaction between the dowel and the timber. 
The foundation modulus calculated in this way can be used 
to predict the stiffness of complete connections for analysis 
of vibration in frames modelled with semi-rigid joints. 
Although the application of the model is limited to vibra¬ 
tion about a non-zero mean load, with no load-sign 
reversal, this form of vibration encompasses various 
important types of in-service vibration of structures, such 
as that induced by turbulent wind or footfall. 

Ein analytisches Modell zur Bestimmung der 
Lochleibungssteifigkeit einer stiftformigen 
Verbindung bei zyklischer Belastung 

Zusammenfassung Bei stiftformigen Verbindungsmit- 
teln ist die Bettung eines Stahldiibels im umliegenden Holz 
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von grundlegender Bedeutung. Die Steifigkeit des Holzes 
im Bereich der Lochleibung wird anhand des Bet- 
tungsmoduls ausgedriickt. Zur Beurteilung der bei zykli¬ 
scher Schwellbelastung auftretenden Sekantensteifigkeit 
wurde ein Standardpriifverfahren zur Bestimmung des 
Bettungsmoduls bei statischer Belastung modifiziert. Es 
wurde gezeigt, dass die stationare Sekantensteifigkeit sig- 
nifikant hoher ist als die statische Steifigkeit bei Anfangs- 
belastung und dass bei einer ausreichend kleinen 
Amplitude der zyklischen Belastungskomponente ein 
einfaches analytisches elastisches Modell zur Bestimmung 
des Bettungsmoduls gut geeignet ist. Das analytische 
Modell basiert auf einer komplexen Spannungsfunktion 
des Holzes im Bereich der Lochleibung sowie der Rei- 
bung zwischen Diibel und Holz. Der so berechnete Bet- 
tungsmodul kann bei den Schwingungsberechnungen 
in Tragwerken mit modellierten nachgiebigen Verbindun- 
gen zur Bestimmung der Steifigkeit dieser Verbindungen 
verwendet werden. Obwohl das Modell nur fur eine 
Schwellbelastung gilt, deckt diese doch verschiedene 
wichtige Schwingungsbeanspruchungen in Bauwerken ab, 
wie zum Beispiel Wind- oder Trittbeanspruchung. 

1 Introduction 

Wind, footfall and machinery can cause vibration in struc¬ 
tures. Vibration above a certain level may be perceived by 
occupants or users of the structure, and this can lead to 
complaint. While such vibration may not necessarily cause 
structural damage, complaint by users represents a ser¬ 
viceability failure, which may result in additional costs in 
modifying the structure to perform appropriately, and may 
also reduce public confidence in a particular form of con¬ 
struction. Such phenomena are magnified by the light 
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weight and flexibility of timber and engineered timber 
products, when compared to alternative structural materials. 

In order to reduce the likelihood of excessive vibration, 
the structural engineer must have accurate design param¬ 
eters. One important parameter is the stiffness exhibited by 
the system. To predict the behaviour of a structure, the 
stiffness properties of each component of the structure 
must be known, and the stiffness of connections can be an 
important factor in the design of timber structures. The 
foundation modulus is the stiffness, per unit length of 
connector, of the resistance of the timber to the embedment 
of a rigid dowel. It can be used as a fundamental parameter 
in assessing the stiffness of dowel-type connections. 
Dowel-type connections are widely used, and include nails, 
screws, bolts and plain dowels. 

In this paper, the foundation modulus applicable to 
vibrations has been investigated, by means of experimental 
work on the interaction between timber and connector 
under dynamic load. The foundation modulus for a single 
dowel has also been theoretically estimated based on the 
orthotropic elastic properties of the timber and the coeffi¬ 
cient of friction between the steel and the timber. Material 
properties for the theoretical model are obtained through 
tests on specimens extracted from the dynamic test pieces, 
and the predictions are tested against the experimental 
results. 

This secant foundation modulus is based on the relative 
deformation of the timber around the dowel for each cycle 
of applied force. It differs from the static foundation 
modulus, which is based on the total deformation of the 
timber under a gradually applied force. The secant foun¬ 
dation modulus is higher than the static value, since the 
latter includes the effect of a low initial stiffness, resulting 
from any lack of fit between the dowel and the hole in the 
timber, and the imperfection of the contact surface between 
dowel and timber. 

This analytical and experimental work focusses on 
vibration without load sign reversal, so, throughout the 
cycle, a compressive force is maintained on the specimen. 
This is an important form of in-service vibration since it 
corresponds, for example, to footfall-induced vibration, 
where the self-weight and imposed load on the structure 
provides a mean force about which the footfall force 
induces oscillation. A similar form of vibration occurs as a 
result of turbulent wind load on a structure, in which the 
mean wind force is greater than the oscillating component 
caused by turbulence. 

2 Background 

As structural engineers strive to create more efficient 
structures, reducing the self-weight of floors and bridges as 


far as possible, serviceability criteria such as vibration 
become an important factor in design. There are examples 
of modern structures in every main construction material 
which have exhibited unacceptable vibration, and timber is 
no exception. The lateral vibration caused by footfall in the 
Lardal footbridge is a notable one (Ingolfsson et al. 2012). 
Heiduschke (2008) noted that the flexibility of connections 
in multi-storey timber buildings means that they may 
experience unacceptable amplitude of vibration under wind 
load, even though they may not reach the height at which 
buildings in conventional structural materials experience 
such problems. 

Dowel-type connections are widely used in timber 
structures, and their relative flexibility means that they play 
an important role in the behaviour of the structure as a 
whole, under both static and dynamic load. Researchers 
have investigated the contribution of dowel-type connec¬ 
tions to the behaviour of trusses (Larsen and Jensen 2000), 
sheathed wall panels (Silih et al. 2005) and frames with 
moment connections (Nishiyama and Ando 2003). Larsen 
and Jensen (2000) identify the conditions under which a 
semi-rigid joint analysis is possible, including the tight fit 
of connectors to avoid unpredictable initial slip. 

Experimental campaigns have been carried out to 
identify the force-displacement response of dowel-type 
connections between timber elements under cyclic loads 
(Allotey and Foschi 2004; Foliente 1995; Foschi 2000; 
Heiduschke 2006a; Noguchi and Komatsu 2004), and the 
results have been used to develop models of the joints. 
These focus predominantly on the post-yield behaviour 
applicable to seismic and extreme wind events. 

The joint model can then be used to form part of a 
numerical model of a complete structure. This type of 
analysis can be used to assess the behaviour of existing 
forms of construction under dynamic loading, as was done 
by He (2001) and Heiduschke (2006b), or to predict the 
performance of proposed devices and techniques, as by 
Loo et al. (2012) for slip-friction connectors, or by 
Awaludin et al. (2008) for pretensioned bolts. 

Small-amplitude vibration of timber structures in service 
differs in two ways from seismic loading: one is that nei¬ 
ther the timber nor the steel yields, though the behaviour of 
the timber cannot be expected to be completely linear. The 
other is that the vibrations are often one-sided. That is to 
say, the vibrations are often small-amplitude oscillations 
around a non-zero mean, with no reversal of the load. This 
is the case for footfall-induced vibration of floors and 
bridges, as well as along-wind vibration of structures due to 
turbulent wind load. Under one-sided vibration, the secant 
stiffness of the connections is independent of the initial 
slip, which has already occurred, and so the condition 
given by Larsen and Jensen (2000) for analysis with semi¬ 
rigid joints is met. 
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Tests on completed timber buildings to investigate their 
behaviour in response to small vibrations are rare. One 
example, for the widely used sheathed stud wall form of 
construction, is Ellis and Bougard’s investigation of the 
BRE Timber Frame 2000 building (Ellis and Bougard 
2001 ). 

This paper seeks to provide a basis for predicting the 
response of timber structures with dowel-type connections 
to small vibrations, based on the fundamental material 
properties and geometry of the connection. The pre-yield 
behaviour of a dowel-type connector in timber is investi¬ 
gated by considering it as a pin-loaded orthotropic elastic 
plate. A complex stress function for an orthotropic plate in 
plane stress, with a hole loaded along its edge, was pro¬ 
posed by Lekhnitskii (1968). This general form of the 
function was adapted to the particular case in which the 
load on the edge of the hole is applied by a rigid circular 
section, referred to as a pin-loaded plate, by De Jong 
(1977). Different forms of the function have been devel¬ 
oped to investigate various phenomena affecting the 
behaviour of the pin-loaded plate. Zhang and Ueng (1984) 
incorporated the effect of friction at the interface between 
pin and plate and Hyer and Klang (1985) allowed for 
deformation of the pin and clearance between the hole and 
the pin, as well as friction. 

Despite the availability of nonlinear finite element 
modelling software capable of dealing with this geometry 
and loading, there continues to be research interest in the 
stress function approach, such as the development of a new 
form of solution by Aluko and Whitworth (2008). One 
reason for the enduring interest in this approach is the 
speed of calculation, and Derdas and Kostopoulos (2011) 
use the stress function to define the loading along the hole 
edge, and apply that loading to a finite element model, thus 
reducing the finite element model to a linear elastic one and 
avoiding the computational cost of carrying out a contact- 
element analysis. 

The potential for application of this method to timber 
has been clear since Lekhnitskii (1968) used the material 
properties for plywood in the example calculations in his 
book. The further research described above was, however, 
generally motivated by the use of pin-jointed composites in 
aircraft and it was Echavarrfa et al. (2007) who noted its 
potential for application to dowel-type connections in 
timber structures, and developed a particular form of the 
stress function for that purpose. 

All of the above researchers used the complex stress 
function approach as a tool for analysis of the stress dis¬ 
tribution around the pin. Since the stress function defines 
the strain field in the plate, it can also be used to consider 
deformation of the connection, and it is on that aspect that 
this research focusses. 


Sjodin et al. (2006) measured the strain field in the 
timber around a group of dowels, and showed that an 
elastic orthotropic plane-stress model for the timber could 
accurately represent the strain distribution. The same 
researchers then showed the influence of friction between 
the dowel and the timber on that strain distribution (Sjodin 
et al. 2008). The interaction between dowel and timber is 
nonlinear not just because of the friction between the two, 
however, but also because the timber surface is initially 
uneven, and this unevenness is compressed under load, 
improving the contact surface. This effect is discussed by 
Dorn et al. (2013), who suggested that the low initial 
stiffness may be due to the compression of the unevenness 
in the contact surface, and noted that under unloading and 
reloading the embedment behaviour is closer to linear 
elastic, as the surface of the timber remains compressed to 
its previous level. In his PhD thesis, Dorn (2012) modelled 
the load-displacement behaviour of the contact surface 
using the overclosure in a contact element analysis of a 
connection, and successfully represented the loading, 
unloading and reloading response. 

This research presents and tests a model for embedment 
in the form of an analytical set of equations which takes 
into account the interface friction and the orthotropy 
elasticity of the timber. The effect of unevenness of the 
contact surface is not taken into account in this model, 
meaning that it is appropriate only to the secant stiffness 
under cyclic loading, in which the contact surface is 
assumed to have been compressed to provide a relatively 
stiff transfer of load to the timber. The model is therefore 
verified against cyclic loading tests. 

3 Materials and methods 

3.1 Analytical method 

The load-deflection response of a circular dowel embed¬ 
ding into a piece of timber depends on the orthotropic 
elastic behaviour of the timber, and also on the frictional 
interaction between the dowel and the timber. The defor¬ 
mation of the cross-section of the steel dowel is assumed to 
be negligible, but the friction between the steel and timber 
determines the distribution of the forces around the face of 
the timber, so must be modelled. 

This method has been developed on the assumption that, 
under small-amplitude vibration significantly below the 
embedment strength, the behaviour of the timber can be 
treated as linear, and therefore modelled using an ortho¬ 
tropic linear elastic model. The appropriateness of this 
assumption is discussed later in this paper, with reference 
to the experimental results. 
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In the model, the foundation modulus for the timber is 
estimated using a stress function of the form first proposed 
by Lekhnitskii (1968), and, in particular, as derived by 
Zhang and Ueng (1984). The definition of the stress 
function exactly follows the work of these researchers, and 
the formulae are collated and reproduced in this paper for 
convenience. The work of these researchers is then 
extended to assess the deformations in a specimen by 
superposition of two semi-infinite plate solutions. 

The stress function represents the stresses in a semi¬ 
infinite orthotropic plate with a half-hole in the surface. 
The forces applied by the dowel along the surface of the 
half-hole are defined on the basis of the geometry and the 
coefficient of friction between the steel and the timber. The 
geometry and notation are shown in Fig. 1. 

3.1.1 General form of stress function 

Lekhnitskii (1968) showed that the general form of the 
stress functions for an infinite orthotropic plate with a hole 
is as given by (1) to (5). The solution of a particular 
problem relies on finding the coefficients a n and b n which 
correspond to the distribution of the load on the edge of the 
hole. C ;,2 are transformed coordinates describing the point 
on the plate under consideration. Given the complex stress 
functions 0j and 0 2 , the displacements u in the x direction 
and v in the y direction are given by (6) and (7). i is the 
imaginary unit. 

0 1 and 0 2 are stress functions defined by their rela¬ 
tionship to the direct and shear stresses in the plate: the 
direct stress in the x direction, for example, is given by 
(j x = 2 Re(p\0\ + p\0' 2 ) where 0\ is the derivative of 0 } 
with respect to zj, and 0\ the derivative of 0 2 with respect 
to z 2 - Re() denotes the real part of what is, in general, a 
complex number in the brackets, and pi and n 2 are com¬ 
plex constants derived from the material properties of the 
plate, as in (8) and (9). 



Fig. 1 Geometry for analysis 
Abb. 1 Geometrische Grundlagen 


This work is not concerned with the stresses in the plate, 
but rather the relative displacements. The displacements in 
the plate are given by the integral of the strains, and are 
therefore proportional to the integral of the stresses. The 
displacements are therefore calculated using the undiffer¬ 
entiated stress functions 0j and 0 2 . A full derivation of the 
general form of the stress functions and their relationship to 
stresses, strains and displacements is given by Lekhnitskii 
(1968). 



z\ = x + p x y zi =x+ ti 2 y 
u = 2Re(/? 1 4> 1 + P20f) + U 
v = 2Re(gi#i + q 2 02 ) 


(1) 

( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 
( 7 ) 


u is the displacement relative to a particular fixed point, 
and includes one non-zero constant of integration, U. /iy, 
p 2 , Ph P21 C 1 i and #2 are derived from the material 
properties of the plate material, the timber, as in (8)— ( 13), 
where E 1 is the elastic modulus of the plate material in the 
x direction, E 2 the elastic modulus in the y direction, G is 
the shear modulus, and Vj the Poisson’s ratio. 
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Lekhnitskii (1968) presents a method to allow the stress 
functions to be defined on the basis of the trigonometric 
series defining the displacements if and v' on the hole edge. 
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The general form of those series is as in (14) and (15). 


(7 


iO 


a m and are coefficients determined by 


this series, the standard form of stress function defined 
by (14) and (15), the stress functions for this particular 
the boundary conditions and * denotes the complex case are defined as in (32) and (33). It can be seen that 
conjugate. the friction coefficient p affects the slip c, since it 

appears in the equations for the constants A 1 and B } , 
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Once the coefficients are determined, they can be used to 
generate the stress functions using (16) to (19), where 
D = P\qi ~ PiQ i- Note that cc 0 and /} 0 do not affect the 
stress distribution in the plate. The overall movement of the 
hole, caused by a net force P on the hole, is allowed for by 
A and B. 
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3.1.2 Boundary conditions 


Zhang and Ueng (1984) used the displacement boundary 
conditions described in (20) to (22), for the geometry 
shown in Fig. 1, to define the stress functions. In doing so, 
they assumed that the dowel remained rigidly circular, so 
that the radial displacement of the timber equals that of the 
face of the dowel, and that some slip occurs in the tan¬ 
gential direction, at all points except 0 — 0. 
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u 0 is, therefore, the displacement of the dowel in the 
semi-infinite plate, and c represents the extent to which 
the dowel slips along the surface of the timber. Zhang 
and Ueng (1984) proposed the trigonometric series in 
(23) to (31) to satisfy these boundary conditions. Using 


v 


/ 


c — 1 c+l\ c -\~ 1 

H-I uq sin 20 H—-— uosinAO 


2c 


2c 


B x -A 


1 


Uq 


A 

gP 


1 


(B 1 -Ax) 


71 2Ax (vi — k — nk) — Bx (iq — k — nk) 


A x = (19n+ Unk A- 10 k- 10vi) 

+ //(11 n — 6 nk + 15 k — 15 v 1 ) 

Bx = 10/2(1 — k) + 10/i(3 k — 3vi + 2 nk + n) 


k 



£1 

E 2 


n 



2(k- vi) 


1 - viv 2 k 
8 = —=-+ 


E 2 


G 


<Z>i = AlnCi + 


c - Iqi- ipi c + 1 ipi\ Uo 


2x 


2D 


2D 


c? 


+ 


c + 1 q 2 - ipi uo 


2c 


2D 


cl 


& 2 — Bln^i + 


(c - l)qx ~ ip 1 _ < 2 + 1 ip A uo 

2D c 2D) ^ 


2x 


2 

2 


c + 1 qx ~ ip 1 uq 


2c 


2D f 


(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


Although the general solution by Lekhnitskii (1968) is for a 
hole in an infinite plate, the form of the stress functions 
proposed by Zhang and Ueng (1984) means that the ver¬ 
tical component of the resultant force applied to the plate is 
zero on the y-axis. These functions therefore approximate a 
semi-infinite plate with a half-hole, which is appropriate for 
the geometry in the experimental work described in Sect. 
3.2. 

It should be noted that this is only an approximate 
solution for a semi-infinite plate with a half-hole, because 
the shear stress on the y-axis is not zero. Zhang and Ueng 
(1984) show that this approximation does not affect the 
stress field significantly, and it is therefore assumed that it 
will not significantly affect the vertical deformations. 
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3.1.3 Superposition for displacements 

The deformation of the dowel for a particular applied load, 
and particular support conditions, will, in general, be dif¬ 
ferent from the deformation u 0 . The strain field is inte¬ 
grated to give the deformations in the plate, resulting, for 
this case, in one non-zero constant of integration, U in (6). 
In this paper, the constant is eliminated by superimposing 
the solutions for two holes loaded in opposite directions. 
By placing the second 2h below the first, as shown in 
Fig. 2, a line of zero displacement is created at a vertical 
distance h from the hole. If the holes are sufficiently well 
separated that there is no significant interaction between 
the stresses on the edges of each hole, then this is a rea¬ 
sonable representation of a plate supported rigidly along 
one edge. 

In order to calculate the movement of the dowel, the 
displacement of the timber in the semi-infinite plate is 
calculated in two positions, u 0 at y — 0, x = r and Uj at 
y = 0, x = 2h—r. u 0 can be calculated directly using (26). 
For Uj, the transformed coordinates and ( 2 are calculated 
using (3) and (4). Since the point is on the x-axis, c/ and £ 2 
are equal in both cases. and & 2 are then calculated using 
(32) and (33), for a unit applied force P. uj is then calcu¬ 
lated using (6). 

The displacement of the dowel relative to the fixed 
surface at x — h is then given by ( uj—u 0 )/2 , and the 
predicted stiffness of the specimen can be calculated as 
2/(U]—u 0 ). 



Fig. 2 Superposition of two infinite-plate solutions to model a 

supported edge 

• • 

Abb. 2 Uberlagerung zweier Endlosplattenlosungen zur Modellie- 
rung einen unterstiitzen Randes 


3.2 Experimental method 

The experimental testing consisted of three parts. First, the 
foundation modulus of the timber under dynamic loading 
was found. Then the parameters required for the analytical 
determination of the foundation modulus were investi¬ 
gated: the frictional resistance of the dowel in a half-hole 
and the elastic moduli of the timber in each grain orien¬ 
tation. The timber used was Douglas fir grown in the 
United Kingdom, which was stored for 12 weeks after 
delivery in a controlled environment, with the temperature 
maintained at 18-22 °C, and relative humidity at 60-65 %. 

3.2.1 Foundation modulus 

The experimental method is based on the static test method 
described in the ASTM standard (ASTM D 5764:1997). 
This method uses a piece of timber with a half-hole in its 
surface. In contrast, methods such as the European Norm 
(EN 383:2007) use a piece with a complete hole through 
the centre. Using a timber piece with a half-hole enables 
the dowel to be held rigidly straight, which means that the 
displacement of the dowel is uniform along its length, and 
the displacement of the loading head is an accurate mea¬ 
sure of the deformation of the timber. This is particularly 
important given the high stiffness of the timber under 
dynamic loading, and the correspondingly small deforma¬ 
tions that need to be measured in these dynamic tests. 

A linear variable differential transformer (LVDT) was 
used to measure the displacement of the loading head, 
except when its 2 mm range was exceeded; in which case, 
the built-in displacement transducer in the servo-hydraulic 
loading machine was used. Figure 3 schematically shows 
the test setup and specimen dimensions for parallel- and 
perpendicular-to-grain specimens. 

3.2.2 Static tests 

In order to characterize the forces applied to the dynamic 
test specimens, their embedment strength was predicted 
according to Eurocode 5 (EN 1995), based on their density. 
The average density of 27 of the specimens was measured 
as 485 kg/nr, and the average moisture content of three 
preliminary specimens was obtained by the oven dry 
method according to the European Norm (EN 13183:2002), 
and found to be 14.9 %. On this basis, the adjusted density 
for the reference value of 12 % moisture content, used in 
Eurocode 5, was calculated as 453 kg/m\ After the 
dynamic tests had been carried out, eight of the dynamic 
test specimens were used for moisture content evalua¬ 
tion. The average moisture content was 14.6 %, with a 
coefficient of variation of 0.065. The preliminary mois¬ 
ture content measurement was, therefore, considered 
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Fig. 3 Schematic diagram of test setup with perpendicular-to-grain ( left ) and parallel-to-grain (right) specimens (dimensions in mm) 

Abb. 3 Schematische Darstellung des Versuchsaufbaus mit Priifung rechtwinklig zur Faserrichtung (links) und in Faserrichtung (rechts) 
(Abmessungen in mm) 


Table 1 Embedment strength of specimens 
Tab. 1 Lochleibungsfestigkeit der Priifkorper 


Grain orientation 

Calculated strength 

Mean measured strength 


(N/mnf) 

(N/mm 2 ) 

Parallel 

29.8 

28.3 

Perpendicular 

18.0 

14.8 


sufficiently accurate for estimation of load capacity. The 
calculated embedment strength is shown in Table 1. Four 
static tests were also carried out, two in each grain orien¬ 
tation, to verify the calculated value. The average of those 
results is also shown in Table 1. 

Since it would be the calculated value that would be 
used to design a connection in practice, the calculated 
embedment strength was used as the reference value for the 
forces applied in the dynamic tests. 

3.2.3 Dynamic tests 

Each specimen was subjected to 1,000 sinusoidal cycles of 
force at 1 Hz. The peak force was either 20, 40 or 80 % of 
the calculated embedment strength of the specimen. 

The amplitude of the sinusoid was given by the R- ratio, 
as defined by (34), where F peak is the maximum com¬ 
pressive force in the cycle, and F trough is the minimum 
compressive force. Nominal /^-ratios of 1.2 and 10 were 
tested. 

R = (34) 

Ftrough 

Each combination of grain orientation, peak force and 
7?-ratio was tested, resulting in a total of 12 tests. There 
were three repetitions of each, giving 36 specimens in total. 
The steady-state values of secant stiffness were assessed 
for comparison with the results of the analytical model. 


The secant stiffness is taken as the slope on the force- 
displacement diagram between the points of maximum and 
minimum displacement in each 1 s period. 

3.2.4 Friction tests 

Tests were carried out to estimate the coefficient of friction 
to use in the analytical model. Figure 4 shows the test setup 
in which a vertical force was applied to two specimens, 
with a long dowel between. The vertical load was applied 
by a servo-hydraulic loading machine, which allowed the 
load to be held constant despite any vertical settlement. 
The dowel was pushed horizontally using a hydraulic hand 
jack, and the force required to do so was continuously 
recorded through a load cell. 

It should be noted that the friction force in this test setup 
is in a different orientation to the friction force developed 
during embedment of the dowel into the timber. Tests by 
McKenzie and Karpovich (1968) showed that the grain 
orientation had no significant effect on the magnitude of 



Fig. 4 Friction test on dynamic test specimens 

Abb. 4 Reibungsversuch an dynamisch beanspruchten Priifkorpern 
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Table 2 Dynamic test specimens from which elastic modulus spec¬ 
imens were extracted 

Tab. 2 Dynamisch beanspruchte Priifkdrper, aus denen Prufkorper 
zur Bestimmung des Elastizitatsmoduls entnommen wurden 


Test 

specimens 

Grain 

orientation 

Peak force 

E peak 

% of static 
failure load 

/?-ratio 

F ,/ 

F trough 

Number 

of 

specimens 

1 

Parallel 

20 

1.2 

2 

2 

Parallel 

20 

10 

2 

3 

Parallel 

40 

1.2 

3 

4 

Parallel 

40 

10 

2 

5 

Parallel 

80 

1.2 

3 

6 

Parallel 

80 

10 

1 

7 

Perpendicular 

20 

1.2 

3 

8 

Perpendicular 

20 

10 

2 

9 

Perpendicular 

40 

1.2 

2 

10 

Perpendicular 

40 

10 

2 

11 

Perpendicular 

80 

1.2 

2 

12 

Perpendicular 

80 

10 

2 


either the static or steady sliding friction coefficients. It is 
considered reasonable, therefore, to assume that the fric¬ 
tional behaviour will be similar in both orientations, so that 
this test provides a suitable estimate of the coefficient of 
friction for use in the analytical model. 

Three tests were carried out with the load parallel to the 
grain, and three with the load perpendicular to the grain, 
using pairs of specimens with the same dimensions as those 
used for dynamic tests. 

3.2.5 Elastic moduli 

The principal elastic moduli were determined from tests 
on specimens extracted from some of the dynamic test 
specimens, following the method given in the European 
Norm (EN 408:2010). The specimens for which elastic 
modulus data were obtained are summarized in Table 2. 
Two elastic modulus specimens were taken from each 
dynamic test specimen, one in each grain direction, so 26 
elastic moduli were obtained in each grain direction, 52 in 
total. 

4 Results 

4.1 Experimental results 

4.1.1 Secant stiffness 

The measured stiffness was not constant throughout each 
test. Figure 5a shows the variation of stiffness through the 


test for a typical specimen, in which each point represents 
the measured secant stiffness for one of the 1,000 cycles of 
force applied. The stiffness was observed to increase, with 
the rate of increase reducing throughout the test. For the 
test shown in Fig. 5a, the stiffness had reached a fairly 
constant value beyond approximately the 700th cycle, and 
this was typical. 

Some parallel-to-grain specimens tested at R — 1.2, 
however, still exhibited a significant rise in stiffness at the 
end of the test, and the most extreme case is shown in 
Fig. 5b. The specimen therefore failed to reach a steady- 
state value of stiffness during the test. This had conse¬ 
quences for the comparison of these experimental results 
with analytical predictions, as discussed in Sect. 4.3. 

The surface of the timber in contact with the steel will 
not have been perfectly smooth and tight fitting due to the 
drilling process and the heterogeneous microstructure of 
the timber. Part of this unevenness and lack of fit may have 
been gradually taken up under the applied force, so con¬ 
tinually increasing the stiffness towards the value of the 
tight-fitting case. 

The scatter in the calculated values of stiffness is 
thought to be due to a combination of two effects. One is 
that the finite sampling rate meant that the highest and 
lowest values of force and displacement in each cycle were 
not necessarily sampled. The other is that electrical noise in 
the force and displacement measurements would add a 
random component to each measurement, and, therefore, to 
the calculated stiffness. The scatter is, therefore, thought to 
be due to the measurement process and not indicative of the 
real behaviour of the specimen. 

4.1.2 Effect of amplitude 

Figure 6 shows the comparison of two force-displacement 
diagrams recorded for individual cycles of force on the 
same specimen. In this case, an oscillating force at R = 10 
was applied for 200 cycles after 1,000 cycles at R = 1.2 
had been completed, so each graph is for the same speci¬ 
men. It can be seen that the secant stiffness is lower for the 
higher-amplitude cycle at R = 10. This reduction in stiff¬ 
ness is caused by nonlinear behaviour under low forces. In 
fact, the force-displacement diagram for R = 10 shows 
significant nonlinear behaviour and hysteresis, while the 
behaviour at R = 1.2 is near linear. 

It is thought that unevenness in the contact surface 
between timber and steel or lack of fit between the dowel 
and the hole may be taken up by the pressure on this 
surface, causing the behaviour to become more linear 
with higher forces. When the force is removed, in the 
R = 10 case, the effect of unevenness and lack of fit may 
be seen again at the lower loads, causing a reduction in 
stiffness. 
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Fig. 5 Variation of stiffness for a a peak force of 40 % of the embedment strength, R = 10, perpendicular to grain and b a peak force of 40 % of 
the embedment strength, R = 1.2, parallel to grain 

Abb. 5 Streuung der Steifigkeit a bei einer Spitzenbeanspruchung von 40 % der Lochleibungsfestigkeit, R = 10, rechtwinklig zur Faserrichtung 
und b bei einer Spitzenbeanspruchung von 40 % der Lochleibungsfestigkeit, R = 1,2, in Faserrichtung 


Fig. 6 Secant stiffness for a a 
peak force of 40 % of the yield 
force, R = 1.2 and b a peak 
force of 40 % of the yield force, 
R = 10 

Abb. 6 Sekantensteifigkeit a) 
bei einer Spitzenbeanspruchung 
von 40 % der FlieBgrenze, 

R = 1,2 und b) bei einer 
Spitzenbeanspruchung von 
40 % der FlieBgrenze, R = 10 




4.1.3 Collated results 

In order to obtain a steady-state value of stiffness for 
comparison between the tests, the mean of the stiffness for 
each cycle after the 700th was used. These stiffnesses are 
plotted against the peak applied force in Fig. 7. The figures 
show that the peak applied force varied between repetitions 
of the same test, particularly in tests with a high peak load 
or a high amplitude of the load. This was because of the 
limited gain in the servo-hydraulic loading machine and the 
variability in specimen stiffness. 

Where R = 10 in both the parallel- and perpendicular- 
to-grain tests, it can be seen that the stiffness is generally 
lower than the corresponding tests with R = 1.2. This is 
due to the effect of nonlinearity at low force, as described 
in Sect. 4.1.2. 

There is a large scatter of stiffnesses in all the tests, but 
particularly in the perpendicular-to-grain tests. This was 
assumed to reflect the variation of material properties in the 


timber, and this assumption was tested by comparing these 
measured values of stiffness with predictions based on the 
material properties of the timber in the specimens and the 
coefficient of friction between the dowel and the timber, as 
described in Sect. 4.2. 

In each of the three tests parallel to the grain, with a peak 
force of 80 % of the embedment strength and R = 1.2, the 
displacement sensor was not correctly positioned for the first 
1,000 cycles of force applied, so the results shown are for a 
second set of 1,000 cycles. 

4.1.4 Friction tests 

Six friction tests were carried out: three with the load 
applied parallel to the grain, and three with the load applied 
perpendicular to the grain. Figure 8 shows the plot of 
applied lateral force against lateral displacement for one of 
the friction tests. This highlights the influence of the rate of 
movement of the dowel. In part A of Fig. 8, the lever of the 


Spring er 






































618 


Eur. J. Wood Prod. (2013) 71:609-622 


Fig. 7 Mean secant stiffness 
for a perpendicular to grain tests 
and b parallel to grain tests 
Abb. 7 Mittlere 
Sekantensteifigkeit a) bei 
Versuchen rechtwinklig zur 
Faserrichtung und b) bei 
Versuchen in Faserrichtung 
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Fig. 8 Friction test results—specimen loaded perpendicular to the 
grain at 8.6 kN, 40 % of the predicted yield force 
Abb. 8 Ergebnisse der Reibungsversuche—Priifkorper wurde mit 
8,6 kN, 40 % der angenommenen FlieBgrenze, rechtwinklig zur 
Faserrichtung belastet 

hand jack was pushed down slowly, which meant that the 
applied force dropped away before the constant rate of 
displacement, required for kinetic friction, could develop. 
The dowel therefore moved a short distance each time and 
then stopped, a phenomenon known as stick-slip 
oscillation. 

In part B of Fig. 8, the lever of the hand jack was moved 
sufficiently quickly that kinetic friction could develop, 
after overcoming the initial, higher force of static friction. 
Since there is continuous movement in the dynamic tests, it 
is considered that the coefficient of kinetic friction is 
appropriate as the input to the analytical model. The 
coefficient of friction to be used in the model was therefore 
estimated based on the kinetic friction in each case. 

In some of the parallel-to-grain friction tests, since a 
higher vertical load was applied, it was not possible to 
depress the lever of the hydraulic jack quickly enough to 
achieve the kinetic friction behaviour as in part B. In these 


cases, the lower value of force in the stick-slip oscillation, 
as in part A, was used to estimate the coefficient of friction. 
Figure 8 shows that the lowest force measured in part A is 
a reasonable estimate of the coefficient of kinetic friction. 

The mean ratio of the horizontal kinetic force to the 
vertical force for all six tests was 0.41, with a coefficient of 
variation of 0.1, giving an estimate of the coefficient of 
friction, for use in the model, of 0.21. There was no clear 
effect of grain orientation on the ratio of normal force to 
friction force. 

4.1.5 Elastic modulus 

The mean measured elastic modulus was 10,569 N/mm in 
the parallel-to-grain direction, with a coefficient of varia- 
tion of 0.21, and 322 N/mm in the perpendicular-to-grain 
direction, with a coefficient of variation of 0.38. 

4.1.6 Estimating shear modulus 

The measured values of the principal elastic moduli pro¬ 
vided two of the four parameters required for the analytical 
model. In addition, the Poisson’s ratio and shear modulus 
are required. The Wood Handbook (Forest Products Lab¬ 
oratory 2010) provides mechanical properties for various 
species of timber, including Douglas fir. There are nine 
independent constants: three elastic moduli corresponding 
to the longitudinal, radial and tangential orientations with 
respect to the grain and ring structure; three shear moduli 
and three Poisson’s ratios. The values for Douglas fir are 
summarized in Table 3. 

The mean elastic moduli from these tests were 
10,569 N/mm 2 in the longitudinal direction and 322 N/ 

j 

mm in the perpendicular direction, which, for most 
specimens, was close to tangential. The ratio of the two 
moduli, was, therefore, quite different from that given in 
the Wood Handbook (Forest Products Laboratory 2010). It 
was not, therefore, considered appropriate to obtain the 
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Table 3 Mechanical properties for Douglas fir, taken from the Wood Handbook (Forest Products Laboratory 2010) 

Tab. 3 Mechanische Eigenschaften von Douglasie, entnommen aus dem Wood Handbook (Forest Products Laboratory 2010) 

E l E t E r G lr G lt G rt v lr v lt v rt 

(N/mm 2 ) (N/mm 2 ) (N/mm 2 ) (N/mm 2 ) (N/m 2 ) (N/mm 2 ) 

14740 737 1002 943 1150 103 0.292 0.449 0.390 


shear modulus as a proportion of the longitudinal elastic 
modulus, using a ratio from published results. Another 
rationale was therefore sought to estimate the shear 
modulus. 

Jones (1975) shows that if the elastic compliances of an 
orthotropic material in plane strain are given by s#, as in 
(35), then the apparent shear compliance, s' 66 , for axes at an 
angle 0 to the principle axes, is given by (36). 




4(511 + 522 — 2 ^ 12 ) sin 20 cos 20 
+ 566(sin 20 — cos 20) 2 



For Douglas fir, as described in the Wood Handbook 
(Forest Products Laboratory 2010), the variation of 
apparent shear modulus with angle in the longitudinal- 
tangential plane can be plotted using (36). The result is 
shown in Fig. 9. 

At 0 — 45°, s' 66 is independent of s 66 , and so the 
apparent shear modulus is independent of the principal 
shear modulus. Using the elastic moduli obtained from the 
tests, it was therefore possible to calculate the apparent 
shear modulus at 0 = 45° for each specimen. The ratio of 



Angle (Radians) 


Fig. 9 Variation of apparent shear modulus with angle from the 
principal direction for Douglas fir 

Abb. 9 Statischer Schubmodul von Douglasie in Abhangigkeit des 
Winkels zwischen Kraft- und Faserrichtung 


the maximum and minimum shear modulus from Fig. 9 
was then assumed to apply to the Douglas fir used in the 
tests, and was used to estimate the principal shear modulus 
for each specimen. 

4.2 Analytical results 

4.2.1 Analytical models 

The two elastic moduli obtained from the tests, the Pois¬ 
son’s ratio for the longitudinal-tangential plane taken from 
the literature and the shear modulus estimated as described 
in the previous section were used to create a set of ana¬ 
lytical models for the dynamic test specimens. Using these 
material properties and applying a unit force, the value of 
the two stress functions could be calculated at each point in 
the specimen. Using the stress functions, the displacement 
of every point in the specimen can be calculated. Those 
displacements are plotted, for illustrative purposes, in 
Fig. 10. 

No free edge effects are allowed for on the sides of the 
specimen. The edge distances are considered to be 



Fig. 10 Vertical displacements in mm in the timber for a typical 
specimen in parallel-to-grain loading (1 kN applied force) 

Abb. 10 Vertikale Verformung im Holz in mm fur einen typischen in 
Faserrichtung belasteten Prtifkorper (1 kN aufgebrachte Last) 
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sufficient that the other free edges will not significantly 
influence the stress or deformation fields. 

4.2.2 Sensitivity to friction coefficient 

The sensitivity of the analytical model to the friction 
coefficient was assessed by varying the coefficient between 
the extremes of values published by McKenzie and Kar¬ 
povich (1968), from their study of friction between steel 
and timber. For comparison with their work, the rate of 
movement between the surfaces in these tests corresponds 
to their classification ‘slow’, and the surface of the steel 
dowel falls between his classifications of ‘rough’ and 
‘smooth’. For all the species of timber tested, the friction 
coefficient in these classifications ranged from 0.08 to 0.62. 
Using the material properties for one of the parallel-to- 
grain dynamic test specimens, the sensitivity analysis is 
summarized in Table 4, comparing these extreme values to 
the mean value obtained from the friction tests and used in 
analysis, 0.21. The resulting variation in stiffness is small 
in comparison to the variability in material properties in the 
timber. It is, therefore, considered that the sensitivity of 

Table 4 Effect of friction coefficient on modelled stiffness 
Tab. 4 Einfluss des Reibungskoeffizienten auf die modellierte 


Steifigkeit 



Friction coefficient 

h 

Stiffness(kN/mm) 

Change from /j. = 0.21 

0.08 

106.2 

3.9 % 

0.21 

110.5 

— 

0.62 

120.8 

9.3 % 


stiffness to friction coefficient is not sufficient to justify a 
more extensive investigation of friction than has already 
been carried out. 

4.2.3 Comparison with experimental results 

The mean and range of analytical results was compared 
with the mean and range of experimental results for each 
grain orientation, and for each R- value. The predicted 
values of stiffness are independent of the magnitude of the 
force applied due to the assumption of elastic behaviour. 
The results for the three values of peak force were, there¬ 
fore, combined for this comparison of predicted and mea¬ 
sured stiffness. Figure 11 plots the range of stiffness and 
the mean stiffness for the experimental results and the 
analytical model, with the results organized by grain ori¬ 
entation and R-v alue. 

4.3 Discussion 

Where the amplitude of vibration is smallest, in the tests 
with R = 1.2, the range of predicted values and measured 
values are similar, as are the mean values for each set. This 
suggests that the elastic model is appropriate for prediction 
of the foundation modulus for small-amplitude vibration. 
For the perpendicular-to-grain tests, the mean predicted 
stiffness is within 5 % of the mean measured stiffness. For 
the parallel-to-grain tests with R = 1.2, the measured val¬ 
ues are lower than predicted. The mean measured stiffness 
is approximately 14 % lower than the mean predicted value 
in these tests. As noted in Sect. 4.1.1, for some of these 
tests, a steady-state value of stiffness was not achieved. 


Fig. 11 Collated parallel- and 
perpendicular-to-grain results— 
the range is shown by the 
vertical line and the mean by 
the horizontal line 
Abb. 11 Zusammengefasste 
Ergebnisse der Versuche 
rechtwinklig zur Faserrichtung 
und in Faserrichtung—Die 
vertikale Linie gibt den 
Streubereich an und die 
horizontale Linie den Mittelwert 
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The measured stiffness in each cycle continued to rise until 
the end of the test. It is, therefore, possible that, if the tests 
were continued, values of steady-state stiffness closer to 
the predicted values might be reached. 

For the greatest amplitude of load, where R = 10, the 
predicted stiffness is significantly higher than the measured 
stiffness in the dynamic tests both parallel- and perpen- 
dicular-to-grain. This can be explained by the reduced 
secant stiffness in these tests due to nonlinear behaviour at 
low force, as described in Sect. 4.1.2. 

5 Conclusion 

Comparison of the analytical model and the experimental 
results shows that, for small-amplitude vibrations and once 
a steady-state dynamic response has developed, an elastic 
orthotropic model is appropriate to estimate the dynamic 
stiffness of the specimen. There is considerable scatter in 
the measured stiffnesses between the tests, but this is 
reflected in the measured variation in elastic moduli. The 
measured stiffness is lower than predicted in the parallel- 
to-grain tests, and this is thought to be a result of the 
imperfection of the contact surface between steel and 
timber not being fully compressed under the applied loads. 

The model presents a method to relate the fundamental 
material properties of the timber to its behaviour as part of 
a connection, and this is done through a compact analytical 
equation. The simplicity of the calculation means that it 
may be applied in circumstances in which a more detailed 
numerical model is inappropriate, due the time or compu¬ 
tational expense that would incur. 

When the amplitude of the oscillating force was 
increased, even though the peak force was not increased, 
the specimens behaved in a more nonlinear way. This 
nonlinearity reduced the secant stiffness which was mea¬ 
sured below the predicted elastic stiffness. It would seem 
appropriate to apply some factor in design to allow for the 
fact that, while an elastic analysis is appropriate for small- 
amplitude vibration, a reduced stiffness is expected as the 
amplitude of vibration increases. Definition of such a factor 
would require further tests with a range of amplitudes of 
force. 

The analysis has been applied to a particular geometry 
in this case, in order to match a standard test method for 
determining the foundation modulus for a dowel-type 
connection in timber. In reality, the foundation modulus 
will depend on the geometry of the timber around the 
connector. Further work will seek to apply the analysis 
method to geometries which occur in real structures. 

Acknowledgments Thanks to Will Bazeley and Neil Price for their 
help and guidance in carrying out the experimental work. 
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